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A common practice in many statistical thermodynamics textbooks when 
calculating the partition function, is to replace the relevant statistical sums 
by integrals over a quasicontinuum of levels, right at the outset.^ Since the 
interlevel separation is inversely proportional to the system's size, when doing 
that, they are really assuming an infinite system, and washing away any finite 
size effects. This, of course, simphfies the problem at hand, but at the same 
time obscures the manner in which size effects enter when dealing with finite 
systems. In this Note, we illustrate the role of finite size effects by evaluating 
the main corrections to the canonical partition function of the ideal gas due 
to a large, yet finite, container volume. 

Sufficient conditions under which it is valid to replace sums by integrals 
have been examined by Stutz^ using the Euler summation formula. For the 
case of the one-particle, one- dimensional partition function, Fox"^ rederived 
Stutz results by using the Jacobi transformation. In the spirit of Fox, we 
will express the partition function of a three-dimensional ideal gas enclosed 
in a impenetrable box in terms of complete elliptic integrals of the first kind, 
whose argument obeys a transcendental equation, and will explicitly derive 
the main finite volume corrections to thermodynamical quantities. A dif- 
ferent approach will be followed for the case of an impenetrable spherical 
container, where we will expand the density of states in powers of the en- 
ergy. Surface terms appear in the expansion coefficients, leading to finite size 
effects in the thermodynamics of the system. 
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We will find that, for both, the box and the spherical container, the 
finite volume correction terms to thermodynamical quantities, have the same 
functional dependence on the ratio (thermal wavelength x area/volume). 
Ideal gas in a box: Consider an ideal monoatomic gas of identical 
particles of mass m, enclosed in a box with sides ai, 02 and 03, at temperature 
T. We will assume a gas temperature (density) high (low) enough so that the 
Boltzmann statistics is applicable. In those conditions, the partition function 
is given by^ 

Z = [zf/N\ (1) 



with 



and 



3 00 

^ = n E exp(-a,n,2) (2) 

1=1 ni=l 



3 fX{T)\' 



(3) 



where A(T) = h/\^3mkBT is the de Broglie thermal wavelength and ks, h 
are the Boltzmann and Planck constants, respectively. At room temperature, 
the spacing between two consecutive arguments in the exponential is very 
small (of the order of 10~^° for m = 10^^^ g, ai = 10 cm and T = 300 K). 
However, since ai depends on the inverse of the product temperature x 
{linear size of the system)"^, according to Eq.(|^), there are situations, such 
as small volumes ( "molecular cavities" ) where such spacing, while still small, 
can grow considerably (of the order of 10~^ for T = 1 K and = 10~^ cm). 
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The function we are interested in is the infinite sum 



sia 



IIexp( 



-an 



(4) 



n=l 



which is a special case of the more general type of sum E^=ig" with |g| < 1. 
This sum can be expressed in terms of the theta function^ 



63(0, q) = l + 2g + 2g^ + 2g9 + ... 



From the identity^ {2K/'kY^'^ = 6s{0,q), we have 



(5) 
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where K{k) is the complete elliptic integral of the first kind 
K{k) = /o"/^ da{l - k sin(a)^) and the argument k satisfies the transcen- 
dental equation'' 

K{1 -k)- {a/7v)K{k) = (7) 

Fig.l shows the plot of K{l — k) and {a/TT)K{k) for different values of (cx/tt). 
As soon as ((t/tt) deviates slightly from unity, the intersection of the two 
curves tends to occur very close to either k = or k = 1, which justifies a 
perturbation expansion around those points. 

At high temperatures (or large container dimensions) (ct/tt) << 1 and 
the root of eq.(|^) lies close to A; = 1, which implies, according to eq.(|^), 
K{k) ^ {7i/a)K{0) = 7r^/2cr. After inserting this into eq.(|D, we obtain 



(8) 



This approximation differs slightly from the standard textbook* expression 
s(cr) = (l/2)(7r/cr)^/^. Table I shows a comparison of both expressions with 
the exact result. We see that eq.® shows better agreement up to values of 
cr as high as a ~ 1. After inserting Eq.(|D into Eqs.(|) and (|l|), we have, 



log(Z) 



A^log 



A^log 



1/2 

8 \cr1a2cr3 j N 

1 / 3 \ 1/2 
1 / -n-^ \ 



e 



+ iV^log(l- 

1=1 

3 

-NT. 



TT 



i=l 



(9) 



8 \(TiCr2Cr3 

We are now in position to evaluate the main corrections to the energy 
U = kBT\d\og{Z)/dT)N,v, pressure P = kBT{d\og{Z)/dV)N,T, specific 
heat C = {dU/dT)Ny and entropy S = kelogiZ) + kBT{d\og{Z)/dT)Ny. 
We obtain : 



Energy: 



1 /A(TMy 
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where Uq = (3/2)A^/cbT and A = 2{aia2 + 0103 + 0203) is the area of the box 
and V = 010203 is its volume. 



Pressure : P ~ Pq 
where Pq = NkeT/V. 

Specific Heat: C ~ Cq 
where Cq = {3/2)NkB. 

Entropy: S ^ Sq — 



^ 1 /A(rMy 

V96^ \ V J 
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where 



^0 = {'i/2)NkB + NkB\og 




(14) 



Thus, the main correction due to the finite volume of the container has the 
form (thermal wavelength x area/volume). 

Ideal gas in spherical container: Let us now consider the same ideal 
gas of the preceding section enclosed in a spherical container of radius R. 
The eigenvalue equation is : ^ ji{kR) = with / = 0,1,2,... where k = 
(Sn'^mE/h'^Y^'^ and ji{x) is a spherical Bessel function. For each /, we have 
an infinite number of solutions, indexed by an integer n = 1, 2, ... and for each 
(/, ra) we have an angular momentum degeneracy of {21 + 1). The energies 
are Ein = (/i^/87r^mi?^)X;^, where X/„ is the nth root of ji{x). In the limit 
of large R, we can approximate the X/„ using McMahon's expansion: 

_,_(/"-!) 4(/i -l)(7/i- 31) _ 32(/i-l)(83/i2_ 982/2 + 3779) 
~ 8b 3~(86)3 15 (86)5 

64 (/i - 1) (6949/i3 - 153855/x2 + 1585743/i - 6277237) 

losW ^ ^ 

with 12= {21 + 1)2 and b = {n + {l/2))7v. 

Now, instead of using Eq.(|T5|) to calculate directly the atomic partition 
function, we will use it to evaluate the density of states, which is all we need 
to calculate the thermodynamic properties of the system. For a container 
of "simple" shape it is natural to expect that the main contributions to the 
density of states will have the form: 



p{E)^p,E'/^ + P2 + p3E-'/^ 



(16) 
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where each term has the typical form of a three, two, and one dimensional 
density of states. The relative importance of each term in Eq.(0) will de- 
pend upon the "dimensionality" of the container. Thus, for our spherical 
container where "volume" predominates over "area" and "length" , the main 
contribution comes from the term proportional to i?^/^, with a small (con- 
stant) correction term. In this case, it is convenient to rewrite Eq.(|l6D in 
terms of dimensionless parameters Ri,R2: 

where we have kept only the surface correction term. We proceed as fol- 
lows: Starting from Eq.(p!5D, we evaluate numerically the integrated density 
of states (f){E), defined as the number of states with energies less than E. 
Then, we perform a least-squares fit of the resulting histogram to the smooth 
function 

^(E) = R, f + R, (^^] ' E (18) 

obtaining the optimal values for i?i,i?2. The approximate density of states 
is then given by Eq . ([T7|) . 

The atomic partition function is given by 

fOO 

z= dE p(E)e-^^. (19) 
Jo 

By using Eqs.([T7|),(p^) and the numerical values obtained for _Ri,i?2, we 
can write the main corrections to several thermodynamics quantities: 
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Energy: 



U^Uo 



: + 0.0461 (^) 



(20) 



where Uq = {2>/2)NkBT and A,V are the area and volume of the sphere, 
respectively. 



Pressure : 
where Pq = NksT/V. 

Specific Heat: 
where Co = {3/2)NkB. 



1 + 0.0461 



(^) 



1 + 0.02305 (^) 



(21) 



(22) 



Entropy: 



where 



S^So- 0.0693NkB 



So = {3/2)NkB + NkBlog 



2.804 



'SmR^kBT' 
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(23) 



(24) 



The method outlined above for the sphere can be used for any other 
container of "simple" shape, provided the eigenergies are known either ana- 
lytically or numerically. To its conceptual simplicity we have to oppose its 
rather slow numerical convergence, requiring very many energies to obtain 
coefficients with a given accuracy. 

As in the case of the box, the main finite volume corrections for the 
sphere have the form constant x (thermal wavelength x area/volume). For 
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both, the box and the spherical container, the thermodynamical quantities 
Uo, Pq, Co and So do not depend on the shape of the container. The main 
finite size corrections have the same functional dependence on the ratio (ther- 
mal wavelength X area/volume) for both shapes. 
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Captions List 



Fig.l : Complete elliptic integral of the first kind. The continuous line is 
K{1 — k), while the dashed, dotted and dot-dashed curves are {a/n)K{k) for 
(a/vr) = 0.5, 1 and 2, respectively. 

Table I Particle in a box: Comparison of the exact atomic partition function 
s{a) (second column) with approximation (j^) (third column) and with the 
standard approximation for high temperatures (fourth column). 
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Eoo 
n=l' 



1/2 



1 0.38631860 0.38622693 0.88622693 

0.5 0.75331414 0.75331413 1.25331414 

0.1 2.30249561 2.30249561 2.80024956 

0.01 8.36226926 8.36226926 8.86226926 

0.001 27.5249561 27.5249561 28.0249561 

0.0001 88.1226926 88.1226926 88.6226926 

0.00001 279.749561 279.749561 280.249561 
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